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1. Introduction 



This work is devoted to the study of the initial-boundary value problem for 
the continuity equation 

dtu + div (bu) = 0, where b :]0,T[xn ^R'^ and u ■.]0,T[xn ^R. (1.1) 

In the previous expression, O C M*^ is an open set, T > is a real number and 
K* . div denotes the divergence computed with respect to the space variable only. 

j£j I The analysis of (1.1) in the case when b has low regularity has recently drawn 

Qv ' considerable attention: for an overview of some of the main contributions, we 

f^ . refer to the lecture notes by Ambrosio and Crippa [2]. Here, we only quote the 

tH- I two main breakthroughs due to DiPerna and Lions [15] and to Ambrosio [1], 

^P ■ which deal with the case when div 6 is bounded and b enjoys Sobolev and BV 

(bounded total variation) regularity, respectively. More precisely, in [15] and [1] 
the authors establish existence and uniqueness results for the Cauchy problem 
posed by coupling (1.1) with an initial datum in the case when Q, = R . 
S^ I In the classical framework where b and u are both smooth up to the boundary, 

H ■ the initial-boundary value problem is posed by prescribing 

' dtu + diY (bu) = m]0,T[xn 
u = g on r^ (1-2) 

u = u at t = 0, 

where u and g are bounded smooth functions and r~ is the portion of ]0, r[x 50 
where the characteristics are entering the domain ]0,r[xO. Note, however, 
that if b and u are not sufficiently regular (if, for example, u is only an L°° 
function), then their values on negligible sets are not, a priori, well defined. 
In § 3.2 we provide the distributional formulation of (1.2) by relying on the 
theory of normal traces for weakly differentiable vector fields, see the works 
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by Anzellotti [5] and, more recently, by Chen and Prid [8], Chen, Torres and 
Ziemer [9] and by Anibrosio, Crippa and Manigha [3]. 
Our main positive result reads as follows: 

Theorem 1.1. Let Q, C M"^ be an open set with uniformly Lipschitz boundary. 
Assume that the vector field b satisfies the following hypotheses: 

1. be L°°{]0,T[xn;R'^); 

2. div6GL~(]0,r[xf]); 

3. for every open and bounded set f!.^, '^ ^, b € Ll^^{[0,T[; BV{^^,;M.'^)). 

Then, given u G L°°(J7) and g € L°°{T~), problem (1.2) admits a unique 
distributional solution u G L°°(]0,T[xO). 

Some remarks are here in order. First, we recall that F" is a subset of 
]0, T[xdn and we point out that L°°(F-) denotes the space L~(F-, ^^ ® J^'^'^). 

Second, we refer to the book by Leoni [17, Definition 12.10] for the definition 
of open set with uniformly Lipschitz boundary. In the case when fi is bounded, 
the definition reduces to the classical condition that il has Lipschitz boundary. 
This regularity assumption guarantees that classical results on the traces of 
Sobolev and BV functions apply to the set fi, see again Leoni [17] for an 
extended discussion. 

Third, several works are devoted to the analysis of the initial-boundary value 
problem (1.2). In particular, we refer to Bardos [G] for an extended discussion 
on the case when b enjoys Lipschitz regularity, and to Mischler [18] for the case 
when the continuity equation in (1.2) is the Vlasov equation. Also, we quote 
reference [7], where Boyer establishes uniqueness and existence results for (1.2) 
and investigates space continuity properties of the trace of the solution on 
suitable surfaces. The main assumption in [7] is that b has Sobolev regularity, 
and besides this there are the technical assumptions that d'lvb = and that 
fi is bounded. See also the analysis by Girault and Ridgway Scott [16] for the 
case when b enjoys Sobolev regularity and is tangent to the boundary. Note 
that the extension of Boyer's proof to the case when b has BV regularity is not 
straightforward. 

Our approach is quite different from Boyer's: indeed, the analysis in [7] is 
based on careful estimates on the behavior of b and u close to the boundary and 
involves the introduction of a system of normal and tangent coordinates at 50, 
and the use of a local regularization of the equation. Conversely, as mentioned 
above, in the present work we rely on the theory of normal traces for weakly 
differentiable vector fields. From the point of view of the results we obtain, the 
main novelties of the present work can be summarized as follows. 

• We provide a distributional formulation of problem (1.2) under the 
solely assumptions that b G L°°(]0,T[xQ;M.'^) and div6 is a locally 
finite Radon measure, see Lemma 3.3 and Definition 3.4 in § 3.2. Con- 
versely, in [7] the distributional formulation of (1.2) requires that b 
enjoys Sobolev regularity. 

• We establish well-posedness of (1.2) (see Theorem 1.1) under the as- 
sumptions that b enjoys BV regularity, while in [7] Sobolev regularity 
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is required. Note, however, that the main novelty in Theorem 1.1 is 
the uniqueness part, since existence can be estabhshed under the solely 
hypotheses that b £ L°°(]0,r[xJ7;M"') and divb £ L°°(]0,r[xJl) by 
closely following the same argument as in [7], see [10] for the technical 
details. We point out in passing that, for the Cauchy problem, the ex- 
tension of the uniqueness result from Sobolev to BV regularity is one 
of the main achievement in Ambrosio's paper [J]. Also, this extension 
is crucial in view of the applications to some classes of nonlinear PDEs 
like systems of conservation laws in several space dimensions, see the 
lecture notes by De Lellis [12] and the references therein. 

• We exhibit some counterexamples (see Theorems 1.2 and 1.3 and Corol- 
lary 1.4 below) showing that, regardless the orientation of b at the 
boundary, uniqueness may be violated as soon as b enjoys BV regu- 
larity in every open set i}^ compactly contained in fi, but the regularity 
deteriorates at the boundary of 0,. Also, as the proof of Theorem 1.3 
shows, if BV regularity deteriorates at the domain boundary, it may 
happen that the normal trace of b at di} is identically zero, while the 
normal trace of bu is identically 1, see § 3.2 for the definition of normal 
trace of b and bu. 

• In [7, § 7.1], Boyer establishes a space continuity property for the so- 
lution of (1.2) in directions trasversal to the vector field b under the 
assumption that b enjoys Sobolev regularity. Proposition 3.5 in the 
present work ensures that an analogous property holds under BV reg- 
ularity assumptions. The property we establish is loosely speaking the 
following: assume S^ is a family of surfaces which continuously de- 
pend on the parameter r and assume moreover that the surfaces are all 
transversal to a given direction. Then the normal trace of the vector 
field ub on S^ strongly converges to the normal trace of ub on S^^ as 
r —7- tq. 

Here is our first counterexample. In the statement of Theorem 1.2, Tr b denotes 
the normal trace of b along the outward pointing, unit normal vector to dQ, as 
defined in § 3.2. 

Theorem 1.2. Let Q be the set il :=]0, +cx3[xIR^. Then there is a vector field 
b :]0, l[xJ7 ^IR3 gy^^ ^/jg^ 

i) 6g L°°(]0,1[xJ1;IR3). 
ii) div b = 0; 
iii) for every open and bounded set i7* such that its closure 0.^, C Jl, we 

have b £ L^{[0,1[; BV{n^]R^)); 
iv) Tr6 = -1 on]0,l[xdi}; 
v) the initial-boundary value problem 

dtu + div {bu) = in ]0,l[xn 

u = on]0,l[xdn (1.3) 

u = att = 
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admits infinitely many different solutions. 

Some remarks are here in order. First, since the vector field b is divergence- 
free, then any solution of (1.3) is a solution of the transport equation 

dtu + b-Vu = 

satisfying zero boundary and initial conditions. Second, the proof of Theo- 
rem 1.2 uses an intriguing construction due to Depauw [14]. 

Finally, note that property iv) in the statement of Theorem 1.2 states that 
the vector field b is inward pointing at the boundary dQ. This fact is actually 
crucial for our argument because it allows us to build on Depauw's construction. 

When the vector field is outward pointing, one could heuristically expect 
that the solution would not be affected by the loss of regularity of b at the 
domain boundary. Indeed, in the smooth case the solution is simply "carried 
out" of the domain along the characteristics and, consequently, the behavior of 
the solution inside the domain is not substantially affected by what happens 
close to the boundary. Hence, one would be tempted to guess that, even in the 
non smooth case, when Tr 6 > on the boundary the solution inside the domain 
is not affected by boundary behaviors and uniqueness should hold even when 
the BV regularity of b deteriorates at the boundary. The example discussed in 
the statement of Theorem 1.3 shows that this is actually not the case and that, 
even if b is outward pointing at di}, then uniqueness may be violated as soon 
as the BV regularity deteriorates at the boundary. 

Theorem 1.3. Let Q be the set Q :=]0, +oo[xM^. Then there is a vector field 
b ■.]0,l[xn ^M^ such that 
i) 6eL°°(]0,l[xO;]R3); 
ii) div 6 = 0; 
iii) for every open and bounded set Q^ such that its closure i}^, C fl, we 

have b £ L^{[0,1[; BV{n,;R^)); 
iv) Tr6= 1 on]0,l[xdn; 
v) the initial-boundary value problem 

dfU + div (bu) = in]0,l[xQ , , 

u = att = ^^^ 

admits infinitely many different solutions. 

We make some observations. First, by a trivial modification of the proof one 
can exhibit a vector field b satisfying properties i), ii), iii) and v) above and, 
instead of property iv), Tr6 = on ]0, l[xdQ. Hence, even in the case when b 
is tangent at the domain boundary, uniqueness may be violated as soon as the 
BV regularity deteriorates at the domain boundary. 

Second, the proof of Theorem 1.3 does not use Depauw's example [14]. The 
key point is constructing a non trivial solution of (1.4) such that u{t, x) > for 
a.e. (i, x) g]0, 1[xQ and Tr (bu) < (note that Tr6 > by property iv) in the 
statement of the theorem). Heuristically, such solution "enters" the domain 0,, 
although the characteristics are outward pointing at the boundary. 
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Third, since Tr6 = 1 on ]0, l[xc?0, then in the formulation of the initial- 
boundary value problem (1.4) we do not prescribe the value of the solution 
u at the boundary. In the proof of Theorem 1.3 we exhibit infinitely many 
different solutions of (1.4) and in general different solutions attain different 
values on ]0, \\xdQ.. However, by refining the proof of Theorem 1.3 we obtain 
the following result. 

Corollary 1.4. Let il he the set fi :=]0, +oo[x]R^, then there is a vector field 
b : ]0, l[xO — ;■ M'^ satisfying requirements i), . . . , v) in the statement of Theo- 
rem 1.3 and such that the initial-boundary value problem (1.4) admits infinitely 
many solutions that satisfy Tr (bu) = on ]0, l[xdi}. 

The additional condition Tr (bu) = in the corollary can be heuristically 
interpreted as (a weak version of) u = on ]0, l[xdrt. 

We also point out that, again by a trivial modification of the proof, one can 
exhibit a vector field b satisfying properties i), ii), iii) and v) in the statement 
of Corollary 1.4 and, instead of property iv), Tr 6 = on ]0, l[x(90. Also, for 
any given real constant c, one can actually construct infinitely many solutions 
of (1.4) that satisfy Tr (bu) = c on ]0, l[xdn. 

Outline. The paper is organized as follows. In § 2 we recall some results 
on normal traces of vector fields established in [3]. In § 3 we establish the 
uniqueness part of the proof of Theorem 1.1 and the space continuity property. 
In § 4 we construct the counter-examples that prove Theorems 1.2 and 1.3 and 
Corollary 1.4. 

Notation. 

• =Sf": the n-dimensional Lebesgue measure. 

• J^"*: the m-dimensional Hausdorff measure. 

• I-l\—E: the restriction of the measure // to the measurable set E. 

• 1e- the characteristic function of the set E. 

• Q: an open set in R"' having uniformly Lipschitz continuous boundary. 

• L°°{]0,T[xdn) := L°°{]0,T[xdn,^^ (g) J^'^-^), where we denote with 
® the (tensor) product of two measures. 

• div b: the distributional divergence of the vector field b :]0, T[xQ — )• M*^, 
computed with respect to the x £ Q, variable only. 

• Div-B: the standard "full" distributional divergence of the vector field 
B. In particular, when B :]0,T[x^ — >• M'^"'"^, then Divi? is the diver- 
gence computed with respect to the {t,x) e]0,T[xO variable . 

• V(/3 : the gradient of the smooth function ip :]0,r[xri — t- M'^, computed 
with respect to the a; S fi variable only. 

• Tr (b, S): the normal trace of the vector field b on the surface SCO, 
as defined in [3] (see also § 2 in here). 

• Tr6: the normal trace of the vector field b on ]0,T[xdQ, defined as 
in S 3.2. 
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• ■A4cxd(A): the class of bounded, measure-divergence vector fields, namely 
the functions B G L°°(A;M^) such that the distributional divergence 
Divi? is a locally bounded Radon measure on the open set A C W . 

• |x|: the Euclidian norm of the vector x G M . 

• supp p: the support of the smooth function p : M — )■ R. 

• Bii{0): the ball of radius R > and center at 0. 

2. Normal traces of bounded, measure-divergence vector fields 

We collect in this section some definitions and properties concerning weak 
traces of measure-divergence vector fields. Our presentation follows [3, §3]. 

Given an open set A C M^, we denote by A^oo(A) the family of bounded, 
measure-divergence vector fields, namely the functions B E L°°(A;]R ) such 
that the distributional divergence Divi? is a locally bounded Radon measure 
on A. 

Definition 2.1. Assume that A C M^ is a domain with uniformly Lipschitz 
continuous boundary. Let B € TWooC^)? then the normal trace of B on dA can 
be defined as a distribution by the identity 

{Tr {B, d A), (f) =[ V(p-Bdx+ f ipd{BivB) y^eC^(M^). (2.1) 
Ja Jk 

This definition is consistent with the Gauss-Green formula if the vector field 
B is sufficiently smooth. In this case the distribution is induced by the inte- 
gration of -B • n on SA, where n is the outward pointing, unit normal vector to 
dA. 

Lemma 2.2. The distribution defined above is induced by an L°° function on 
dA, which we can still call Tic{B,dA), with 

||TV(S,9A)||i.o(aA)<||i3||L-(A)- (2.2) 

Moreover, if T, is a Borel set contained in 9Ai n9A2 and ifni = n2 on S, then 

Tr(5, Ml) = Tr(5, M2) Jf^"^ - a.e. on S. (2.3) 

Starting from the identity (2.3), it is possible to introduce the notion of 
normal trace on general bounded, oriented, Lipschitz continuous hypersurfaces 
S C R-^. Indeed, once the orientation of ns is fixed, we can find Ai C M.^ such 
that E C (?Ai and the normal vectors n^ and ni coincide. Then we can define 

Tr~(5, S) := Tr{B, Mi). (2.4) 

Analogously, if A2 ^ M is an open subset such that S C 9A2, and ^2 = —ns, 
we can define 

Tr+(B, S) := -TV(5, M2). (2.5) 

Note that we have the formula 

(DivS)Li;= (Tv+{B,T,)-Tr-{B,T,))je^~^\_T,. (2.6) 



In particular, Tr^ and Tr coincide J^^ ^-a.e. on S if and only if S is negligible 
for the measure Div B. 
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We now go over some space continuity results established in [3, §3]. We first 
recall the definition of a family of graphs. 

Definition 2.3. Let / C R be an open interval. A family of oriented surfaces 
{Sr}rG/ ^ ^'^ is a family of graphs if there are a bounded open set D C ]R^~^ 
and a Lipschitz function / : D — ?• M such that the following holds. There is a 
system of coordinates (xi, • • • , xjy) such that, for every r £ I, 

Sr = { (Xi, . . . , Xn) : f{xi, ..., XN-l) - XN = r} 



and Sr is oriented by the normal (— V/, 1)/y^1 + |V/p. 

We now quote [3, Theorem 3.7]. 

Theorem 2.4. Let B G A^oo(I^^) o^c^ ^ei {S^jre/ ^e a family of graphs as in 
Definition 2.3. Given vq G /, we define the functions ao,ar : -D — )■ M by setting 

ao{xi,... ,XAr_i) := Tv'{B,T.rJ{xi, . . . ,XN-1, fixi, . . . ,XN-l) -To) 

and 

arixi,...,X]\f_i) := Tr + (5,S,.)(xi,...,XAr_i,/(xi,...,XAr-i) -r). 
Then we have 

ar ^ ao weakly* in L°°(D,^^-^LL>) as r ^ r+ . 

3. Proof of Theorem 1.1 

3.1. Preliminary results. In this section we establish some results that are 
preliminary to the distributional formulation of problem (1.2). 

Lemma 3.1. Let B be a locally bounded vector field on M and let {pe}o<e<i 
be a standard family of mollifiers satisfying supp/3^ C ^^(O) for every e 6]0, 1[. 
The divergence of B is a locally finite measure if and only if for any K 
compact in M there exists a positive constant C such that the inequality 

\\D\wB*p,\\li(k)<C (3.1) 

holds uniformly in e g]0, 1[. 

Proof. If Div i? is a locally finite measure the inequality (3.1) is satisfied on any 
compact K for some constant C independent from e. 

On the other hand, the sequence (Divi?) * ps = Div (B * p^) converges to 
Div B in the sense of distributions and the uniform bound (3.1) implies that we 
can extract a subsequence which converges weakly in the sense of measures. D 

Lemma 3.2. Let A C M^ be an open subset with uniformly Lipschitz continu- 
ous boundary and let B belong to A^oo(A). Then the vector field 

' B{z) zgA 



' U otherwise 



belongs to A^oo(I 



iN\ 
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Proof. We only need to check that the distributional divergence of i? is a locally 
bounded Radon measure. Given e g]0, 1[ we define the e- neighborhood of dA 
as 

dA, = {z£R^ : dist(z,M) < e}. 

Any compact subset K of M can be decomposed as follows: 

K ={Kr\{A\ dAe)) u{Kn dAe) u{K\{AU dAe)) . (3.2) 

Also, note that Div {B * p,) is zero on i^ \ (A U SA^) and that its L^ norm is 
uniformly bounded on A' n (A \ dAg,). Moreover, 

/ |div {B * p,)\dz= \B* Vpe\ dz 

JKndA, JKndKe (3.3) 

<PllL-(R'V)||Vpe||Li(RiV)^'^(Kn5Ae). 

We observe that ||i3||2,oomiV) = ||i?||/^oo(^), that J^^ [K r\dA,) < C*e and that 



|L°°(A)) 
II^PellLi{R^) ^ C'**/^ for suitable constants C^, > and C** > 0. Hence, 

/ |div {B * pe)\ dz < ||i?||£^oo(A)C*C** 

JKndAe 

and by relying on Lemma 3.1 we conclude. D 

3.2. Distributional formulation of problem (12). We can now discuss the 
distributional formulation of (1.2). The following result provides a distribu- 
tional formulation of the normal trace of b and bu on ]0,T[xdil.. 

Lemma 3.3. Let fi C M*^ 6e an open set with uniformly Lipschitz boundary and 
let T > 0. Assume that b G L°°{]0,T[xQ;M.'^) is a vector field such that div 6 
is a locally finite Radon measure on ]0,T[xil. Then there is a unique function, 
which in the following we denote by Tib, that belongs to L'^{]0,T[xdQ) and 
satisfies 

/ / TrbipdJ^'^-'^dt- ip{0,x)dx= dttp + b ■ V (p dxdt 

Jo Jan Jn Jo Jn /o a\ 

+ [ [ tpdidivb) y^eC^{[0,T[xM.'^). 
Jo Jn 

Also, ifw€ L°°(]0,r[xO) and f G L~(]0,r[xO) satisfy 

I [ w{dtV + b-V7])dxdt+ [ [ f7]dxdt = yr]eC^i]0,T[xn), (3.5) 
Jo Jn Jo Jn 

then there are two uniquely determined functions, which in the following we 

denote by Tr (bw) G L'^{]0,T[xdn) and wq £ L°°(0), that satisfy 

/ / Tr {bw)ipdJl^'^-^dt- / ip{0,-)wodx 
Jo Jen Jn 

= [ [ w{dt^ + b-V^)dxdt+ [ [ fipdxdt V^? G C~([0,r[xM'^), 
Jo Jn Jo Jn 

(3.6) 
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Note that requirement (3.5) is nothing but the distributional formulation of 
the equation 

dtw + div (bw) = f m]0,T[xn. (3.7) 

Also note that the existence of the function wq follows from Lemma 1.3.3 in 
[11], the new part is the existence of the function Tr (bw). 

Proof. We first establish the existence of a function Tr6 satisfying (3.4). Note 
that the uniqueness of such a function follows from the arbitrariness of the test 
function ip. We define the vector field B : M"'+^ — )• M"'+-'^ by setting 

S(t,x):=| Q elsewhere in M'^+i ^^'^^ 

and we note that DivSl, „r „ = div6, therefore B satisfies the hypotheses of 

Lemma 3.2 provided that A :=]0,r[xJ7. Hence, B G Moo(^'^^'^)- We apply 
Lemma 2.2 and we observe that Tr {B,dA)\f . „ = —1. We can then conclude 
by setting 

TV6:=Tr(i3,aA)|^„^^j^,^. 

The existence of the function Tr (bw) satisfying (3.6) can be established by 
setting 

^ ' elsewhere m ™'^+^ ^ ' 



and observing that condition (3.5) implies that DivCL^r ^ = /. We can then 
conclude by using the same argument as before, by setting 

wo:=-TT{C,dA)\^^^^^ and Tr {bw) := Tr {C, 8 A)]^^.^^^^^. (3.10) 

D 

We now state the rigorous formulation of problem (1.2). 

Definition 3.4. Let fi C M'^ be an open set with uniformly Lipschitz boundary. 
Assume that b E L°°{]0,T[xQ;W^) is a vector field such that divb is a locally 
finite Radon measure on ]0, T[xQ. A distributional solution of (1.2) is a function 

u£ L°°{]0,T[xn) such that 

i) u satisfies equation (3.5) with / = 0; 
ii) Wq = u; 
iii) Tr (bu) = gTib on the set T~ which is defined as follows: 

r- := {{t,x) e]0,T[xdn : (Tr b){t,x) < O}. 

Note that in Definition 3.4 we only assume / = for the sake of simplicity. 
By removing the condition / = from point i) we obtain the distributional 
formulation of the initial-boundary value problem obtained by replacing the 
first line of (1.2) with (3.7). 
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3.3. Proof of Theorem 1.1. First, we observe that the existence of a solution 
of (1.2) is established in [10] by closely following an argument due to Boyer [7]. 
The argument to establish uniqueness is organized in two main steps: in § 3.3.1 
we show that, under the hypotheses of Theorem 1.1, distributional solutions 
of (1.2) enjoy renormalization properties. Next, in § 3.3.2 we conclude by 
relying on a by now standard argument based on Gronwall Lemma. 



3.3.1. Renormalization properties. We fix u distributional solution of (1.2) and 
we proceed according to the following steps. 

Step 1: we use the same argument as in Ambrosio [1] to establish renormal- 
ization properties "inside" the domain. More precisely, the Renormalization 
Theorem [1, Theorem 3.5] implies that the function u^ satisfies 

/ u^{dtip + b-V%l))dxdt- I I u^d\\biljdxdt = - I u^'ip{0,-)dx 
Jn Jo Jn Jn 

yi^ec^{[o,T[xn). 

(3.11) 

Step 2: we establish a trace renormalization property. 

First, we observe that by combining hypothesis 3 in the statement of Theo- 
rem 1.1 with Theorem 3.84 in the book by Ambrosio, Fusco and Pallara [4] we 
obtain that the vector field B defined as in (3.8) satisfies B{t, ■) £ BV{Q^:) for 
every open and bounded set fi* C M'^ and for .=Sf^-a.e. t £]0,r[. 

Next, we recall that the proof of Lemma 3.3 ensures that the vector field uB 
belongs to Aioo{^'^~^^)- We can then apply [3, Theorem 4.2], which implies the 
following trace renormalization property: 



Ti{u^b)it,x) = < 



Tr 6 y 


Tr b{t, x) / 





Tr b{t,x) = 0. 



(3.12) 



Some remarks are here in order. First, to define Tr (u^6) we recall (3.11), use 
Lemma 3.3 and set 

Tt{u^) := Tr (n2i?,M)|]o_^j^g^, (3.13) 

where A=]0,r[xJ]. 

Second, note that, strictly speaking, the statement of [3, Theorem 4.2] re- 
quires that the vector field B has BV regularity with respect to the {t, x)- 
variables, which in our case would imply some control on the time derivative of 
b. However, by examining the proof of [3, Theorem 4.2] and using the particular 
structure of the vector field B one can see that only space regularity is needed 
to establish (3.12). 



IBVP FOR THE CONTINUITY EQUATION 11 

Step 3: by combining (3.11) with (3.13) and recalling Lemma 3.3 we infer 
that 

[ [ TT{uH)ipdM"^-^dt- [ u^ip{0,x)dx= I [ u^ {dt^ + b ■ V if) dxdt 
Jo Jdn Jn Jo Jn 

[ u^divb^dxdt yip£C^{[0,T[xR'^). 
Jn 

(3.14) 

3.3.2. Conclusion of the proof of Theorem 1.1. We conclude by following a by 
now standard argument, see for example the expository work by De Lellis [13, 
Proposition 1.6]. We proceed according to the following steps. 

Step A: we observe that, the equation being linear, establishing that the 
distributional solution of (1.2) is unique amounts to show that, if n = and 
^ = 0, then any distributional solution satisfies u = 0. Also, note that in the 
remaining part of the proof we use [11, Lemma 1.3.3] and we identify u^ with 
its representative satisfying that the map t i— )• u^(t, •) is continuous in L^^{Q.) 
endowed with the weak* topology. 

Step B: we fix t g]0, T[ and we construct a sequence of test functions ipn as 
follows. First, we choose a function /i : [0, +oo[— )• M such that 

h e C^{[0, +oo[), h>Oandh' <0 everywhere in [0, +oo[. (3.15) 

Next, we set 

iy{t,x) :=h{\\b\\Lo^\t-i\ + \x\) (3.16) 

and we observe that v satisfies 

dti^ + b-Viy<dti^+\\b\\Lo^\Viy\<0 ^'^+^-a.e. (t, x) e] - oo, t[xR'^. (3.17) 
We then choose a sequence of cut-off functions Xn G C^([0, +oo[) satisfying 

Xn = 1 on [0, i], Xn = on [t + l/n,+oo[, Xn ^ everywhere on [0, +oo[. 

(3.18) 
Finally, we set 

Mt, x) ■■= XnitHt, x) (t, x) G [0, T[xR'^ 

and we observe that tpn > everywhere on [0,+oo[xIR'^ and that tpn is com- 
pactly supported in [0,T[x]R'^ provided that n is sufficiently large. 

Step C: we use ipn as a test function in (3.14). First, we observe that by 
recalling that 3 = and u = and by using the renormalization property (3.12) 
we obtain that the left hand side of (3.14) is nonnegative, namely 

0< / u ux'ndxdt + / / Xn^ [dtv + b ■\/v)dxdt — I I div6u vxndxdt. 
Jo Jn Jo Jn Jo Jn 

Next, we let n — t- -|-oo and by recalling properties (3.17) and (3.18) we obtain 

i'{i,-)u (t,-)dx< ||div6||ioo / / uu dxdt. 
n Jo Jn 
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We can finally conclude by using Gronwall Lemma and the arbitrariness of the 
function h in (3.16). This concludes the proof of Theorem 1.1. D 

3.4. Rigorous statement and proof of the space continuity property. 

We provide a rigorous formulation of the analogue of the space continuity prop- 
erty stablished in the Sobolev case by Boyer in [7, § 7.1]. 

Proposition 3.5. Let b he as in the statement of Theorem 1.1, u G L°°(]0,T[xil) 
he the distrihutional solution of (1.2) and B G A^oo(I^'^^^) be the same vector 
field as in (3.8). Given a family of graphs {S^lrG/ ^ 1^°' as in Definition 2.3, 
we fix rQ £ I and we define the functions 70, 7r '■]0,T[xD -^M by setting 

7o(t, xi, . . . , Xd-i) := Tr~{uB, ]0, T[xT,ro){t, xi, . . . , Xd-i,f{xi, ..., Xd-i) - ro) 

and 

'yr{t,xi,...,Xd-i) := Tv~^{uB,]0,T[xT,r){t,xi,. . . ,Xd^i, fixi, . . . ,Xd-i) -r). 

Then 

Ir — ^ 7o strongly in L {]0,T[xD) as r ^ r^. (3.19) 

Proof. The argument is organized in three steps. 

Step 1: we make some preliminary considerations and introduce some no- 
tation. With a slight abuse of notation, we consider 6 as a vector field defined 
on M'^+^j set equal to zero out of ]0,T[xi}. 

By combining hypothesis 3 in the statement of Theorem 1.1 with [1, Theorem 
3.84] we obtain that b{t, •) G BViocO^'^) for ^^-a.e. t G R. Hence, the classical 
theory of BV functions (see for instance [4, Section 3.7]) ensures that the outer 
and inner traces b{t, •)^ and b{t, •)^ are well-defined, vector valued functions 
for ^"^-a.e. i G M and for every r. 

Step 2: given B as in (3.8), we define the functions /3o>/3r :]0,r[xL' — ;> R 
by setting 

/3o{t,xi,... ,Xd-i) ■■= Tr~{B,]0,T[xT,ro){t,xi, . . . ,Xd-i, f{xi, . . . ,Xd-i) - ro) 

and 

/3r{t, xi,..., Xd-i) := Tr +(S, ]0, T[x S^) {t, xi, . . . , Xd-i, f{xi, ..., Xd-i) - r) . 

We claim that 

(5r -^ /3o strongly in L^(]0, T[xD) as r -^ r+. (3.20) 

To establish (3.20), we first observe that by using [i, Theorem 3.88] and an 
approximation argument one can show that for every r G / we have 

/3r = b'^^-rh, and l3o = b^^ ■ rh ior ^'^-a.e. {t,x) £]0,T[xD. 



In the previous expression, fh = (— V/, 1)/ y'T+|V/P is the unit normal vector 
defining the orientation of T,r. Also, by again combining [4, Theorem 3.88] with 
an approximation argument we get that 

' \f3r-f3o\dxi...d^,_,dt< [ \Db{t,-)\{S)dt, 

D Jo 
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which imphes (3.20). In the previous expression, \Db{t,-)\ denotes the total 
variation of the distributional derivative of b(t, •), and S is the set 

S := ^{xi,...,Xd-i,Xd) :{xi,...,Xd-i) e D and 

f{xi,...,Xd-i) -r < Xd< f{xi,...,Xd-i) -roj . 

Step 3: we conclude the proof of Proposition 3.5. First, we observe that 
due to Theorem 2.4 we have that 

7r ^ 70 weakly in L'^{]0,T[xD) as r ^ r+. (3.21) 

Next, we recall that -jr is the normal trace of uB and that /3r is the trace of B, 
so that by applying [3, Theorem 4.2] we get 

7^ = /3,.Tr + (u^S,]0,r[xS^) and -y^ = f3oTi^ {u'^B,]0,T[xT.ro). (3.22) 

By combining (3.20) with the uniform bound ||/3r||Loo < ||6||loo we infer 
that f3r —7- /3o strongly in L'^(]0,T[xD). Then we apply Theorem 2.4 to 
Tr ~^ {u"^ B ,]0,T[xT,r) and hence by recalling (3.22) we conclude that 

1r -^ To weakly in L^{]0,T[xD) as r ^ r+. (3.23) 

By using (3.21), we get that (3.23) implies that ^r — ^ 7o strongly in 
L'^{]0,T[xD) and from this we eventually get (3.19). D 

4. Counter-examples 

4.1. Some notation and a preliminary result. For the reader's conve- 
nience, we collect here some notation we use in this section. 

• Throughout all § 4, $7 denotes the set ]0,+oo[xR^. 

• We use the notation (r, y) g]0, -|-oo[x]R^ or, if needed, the notation 
{f^yi-,y2) S ]0,-|-oo[ X M X M to denote points in 0. 

• div denotes the divergence computed with respect to the (r, y)-variable. 

• Div denotes the divergence computed with respect to the (t, r, y)- 
variable. 

• divy denotes the divergence computed with respect to the y variable 
only. 

• We decompose ]0, l[x$7 as ]0, l[x$7 = A+ U A~ U 5, where 

A+ :={(t,r,y) G]0,l[xrj : r >t} (4.1) 

and 

A-:={(t,r,y)G]0,l[xO: r < t}, (4.2) 

while S is the surface 

5:={(i,r,y)G]0,l[xO: r = t}. (4.3) 

We also observe that, thanks to [11, Lemma 1.3.3], up to a redefinition of n(t,3;) 
in a negligible set of times, we can assume that the map 1 1— t- n(t, •) is continuous 
from ]0, 1[ in L°°{Vt) endowed with the weak-* topology, and in particular 

n(t, •) ^ no in L~(f7) as t -^ 0+, 
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where uq the value attained by u at i = 0, as in Lemma 3.3. 

4.2. Proof of Theorem 1.2. The proof is organized in three steps. 

Step 1: we recall an intriguing example due to Depauw [14] which is piv- 
otal to our construction. In [14], Depauw explicitly exhibits a vector field 
c :]0, +oo[x]R^ —7- M? satisfying the following properties: 

a) cG L°°(]0,+oo[xM2;m2). 

b) For every r > 0, c{r, •) is piecewise smooth and, for almost every y G M^, 
the characteristic curve trough y is well defined. 

c) divyC = 0. 

d) c€ Ll^{]0,l[;BVioc{^^;R^)),hntc^ L\[0,l[;BVioc{R^;R^)). Namely, 
the BV regularity deteriorates as r — )• O"*". 

e) The Cauchy problem 

( drw + dwy{cw) = on]0, 1[xM2 

\u; = atr = ^ ' ' 

admits a nontrivial bounded solution, which in the following we denote 
by v{r,y). 

Step 2: we exhibit a vector field b satisfying properties i), . . . , v) in the 
statement of Theorem 1.3. We recall that the sets A"*", A~ and S are defined 
by (4.1), (4.2) and (4.3), respectively. We define the vector field b :]0, l[xQ. — )■ 
M^ by setting 

7 /i N f (l)C(r, u)) in A~ ,, _, 

^(*'^'^)^=| (lio) ^ inA+ ^^-'^ 

In the previous expression, c is Depauw's vector field as in Step 1. By relying 

on properties a), c) and d) in Step 1 one can show that b satisfies properties 

i), ii), iii) in the statement of Theorem 1.2. 

Next, we recall that the initial-boundary value problem (1.3) admits the 
trivial solution u = and that the linear combination of solutions is again 
a solution. Hence, establishing property v) in the statement of Theorem 1.2 
amounts to exhibit a nontrivial solution of (1.3). We define the function u by 
setting 

/, \ f v{f,v) in A~ ,, „■. 

^(*'^'^)^=|o inA+, (^-6) 

where v is the same function as in Step 1. 

Step 3: we show that the function n is a distributional solution of (1.3). We 
set C := {u,bu) and we observe that by construction DivC = on A"*". Also, 
property e) in Step 1 implies that DivC = on A". Finally, by recalling (2.6) 
we infer that DivCL5 = since the normal trace is on both sides. 

We are left to show that the initial and boundary data are attained. First, 
we observe that u{t, •) ^ as t — )• 0"'' and hence uq = by the weak continuity 
of u with respect to time. Next, we fix an open and bounded set D CM? and 
we define the family of graphs {S^lrgicif ^]0) l[xO by setting 

S, :={(t,r, 2/1,2/2) : t G]0, 1[ and (2/1,2/2) G D}. (4.7) 
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The orientation is given by the vector (0, —1,0,0). We point out that require- 
ment e) in Step 1 imphes that v{r, •) ^ as r — t- 0"*". Hence, by recalhng that 
b is given by (4.5), we obtain that Tr"'"(C, S^) ^ as r — )• 0"*". By recahing 
Theorem 2.4 and Definition 3.4, we infer that Ti (bu) = on ]0, l[xdO,. This 
concludes the proof of Theorem 1.2. □ 

4.3. Proof of Theorem 1.3. The proof is divided in four main steps: 

(1) in § 4.3.1 we construct the auxihary vector field /5fc, which will serve as 
a "building block" for the construction of the vector field b; 

(2) in § 4.3.2 we define the vector field b; 

(3) in § 4.3.3 we show that b satisfies properties iii) and iv) in the statement 
of Theorem 1.3; 

(4) finally, in § 4.3.4 we exhibit a non trivial solution of the initial-boundary 
value problem (1.4). Since the problem is linear, any linear combination 
of solutions is also a solution and hence the existence of a nontrivial 
solution implies the existence of infinitely many different solutions. 

4.3.1. Construction of the vector field /3k. We fix /c G N and we construct the 
vector field /3fc, which is defined on the cell 

(r,2/i,y2)G]0,4-2-'^[x]0,4-2-'=[x]0,4-2-'=[. 

We split the r-interval ]0,4 • 2~'^[ into four equal sub-intervals and we proceed 
according to the following steps. 

Step 1: if r g]0,2~'^[, we consider a "three-colors chessboard" in the {yi,y2)- 
variables at scale 2~^ as in Figure 1, left part. The vector field (3^ attains 
the values (1,0,0), (—5,0,0) and (0,0,0) on dashed, black and white squares, 
respectively. Note that [3^ satisfies 

dw(3k = on]0,2-'=[x]0,4-2-^'[x]0,4-2-'=[ (4.8) 

since /3fc is piecewise constant and tangent at its discontinuity surfaces. 
Here is the rigorous definition of Pk'- we set 

Dk:= U ](2n)2-^(2n + l)2-'=[x](2m)2-^(2m + l)2-^[ 

n,m=0,l 

and 

Bk:= \J ](2n + l)2-^(2n + 2)2-^•[x](2m + l)2-^(2?T^ + 2)2-'=[. (4.9) 

n,m=0,l 

Note that D^ and B^ are represented in the left part of Figure 1 by dashed and 
black regions, respectively. Next, we define 

r (1,0,0) if(yi,y2)Gl?fc 



(3kir,yi,y2) := 



(-5,0,0) if{yi,y2)eBk (4.10) 



(0,0,0) elsewhere on ]0, 4 • 2"'^ [ x]0, 4-2-'' [ 
Step 2: if r g]2~'^,2 • 2~^[, then the heuristic idea to define /3fc is that we 
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Figure 1. The vector field /3fc(?^, yi,2/2) for different values of 
r: the dashed, black and white squares are the region where /3fc 
attains the values (1,0,0), (—5,0,0) and (0,0,0), respectively. 



2/2 



2/2 



2/2 





re]0,2-''[ 



r = 2 • 2- 



2/1 2/1 

r e]3-2-'=,4-2-''^[ 



want to (i) horizontally leftward slide the rightmost dashed squares and (ii) 
horizontally rightward slide the leftmost black squares. The final goal is that 
at r = 2 • 2"^^ we have reached the configuration of the vector field described in 
Figure 1, center part. The nontrivial issue is that we also require that 



div/3fc = on ]0,2 • 2-'=[x]0,4 • 2-*^[x]0,4 • 2-'=[. 



(4.11) 



To achieve (4.11), we employ the construction illustrated in Figure 2: the vector 
field /3fc attains the value (1,0,0) on the horizontal part of the dashed region, the 
value (1, —1, 0) on the inclined part of the dashed region and the value (0, 0, 0) 
elsewhere. Note that (4.11) is satisfied because /3fc is piecewise constant and 
it is tangent at its discontinuity surfaces on the interval r s]2~'^,2 • 2~^[. We 
conclude by recalling (4.8) and by observing that the normal trace is continuous 
at the discontinuity surface r = 2~ and hence no divergence is created there. 
Here is the rigorous definition of /Sj. for r G]2~'^, 2 • 2~^[: 

{ (1,0,0) if(2/i,2/2)G]0,2-'=[x]0,2-^^[ 

or(?/l,?/2)G]0,2-'=[x]2•2-^3•2-^[ 

(1,-1,0) if-r + 3-2-'= <yi < -r + 4-2-'^ 

and 2/2 G]0, 2-''[ or ya G]2 • 2-'=, 3 • 2-*^[ 



/3kir,yi,y2) ■■-- 



(-5, 0,0) if (yi, y2) G]3 • 2-^4 • 2-^[x]2-\ 2 ■ 2-^[ 

or (yi, y2) G]3 • 2-^ 4 • 2-'=[x]3 • 2-^ 4 • 2-^[ 

(-5,-5,0) if r <yi <r + 2-'= 

andya G]2-*^, 2 • 2-'=[ or ya G]3 • 2-'=,4 • 2-'=[ 



(0,0,0) elsewhereon]0,4-2-'=[x]0,4-2-'=[ 



(4.12) 



Step 3: if r g]2 • 2 '^,3-2 ^[, the heuristic idea is defining /3fc in such a way 
that (i) we push up the lower black region in Figure 1, central part, (ii) we 



IBVP FOR THE CONTINUITY EQUATION 



17 



Figure 2. The vector field Pk{r,yi,y2) for y2 G]0,2 ^■[ and 
7/2 e]2 • 2-^=, 3 • 2-'=[. The field /3fc attains the value (1,-1,0) 
in the inclined dashed region. 



yi 



\ 



/• = 2 ■ !'-'■ 



pull down the upper dashed region in Figure 1, central part and (iii) we satisfy 
the requirement that /3fc is divergence- free. This is done by basically using the 
same construction as in Step 2. Note that at r = 3 • 2~ we have reached the 
configuration described in Figure 1, right part. 

Here is the rigorous definition of Pk for r ^]2 ■ 2~ , 3 • 2~ [: 

f (1,0,0) if (2/1,7/2) G]0, 2 •2-^^[x]0,2-^'[ 



h(r,yi,y2) 



(1,0,-1) if7/iG]0,2-2-'^[ 

and -r + 4 • 2^^ < 7/2 < -r + 5 • 2^^" 



(-5, 0, 0) if (yi, ya) G]2 • 2-^ 4 • 2-'=[x]3 • 2-^ 4 • 2-'=[ 

(-5,0,-5) if yi e]2-2-'^',4-2-'=[ 
and r — 2~^ < y2 < f 



(0, 0, 0) elsewhere on ] 0, 4 • 2"^ [ x ] 0, 4 • 2"'= [ 

Step 4: if r g]3 • 2~'^,4 • 2~^[, then we consider the "three colors chessboard" 
in the (yi,y2)-variables at scale 2 • 2"^^ illustrated in Figure 1, right part. The 
vector field /3fc attains the value (1, 0, 0), (—5, 0, 0) and (0, 0, 0) on dashed, black 
and white regions, respectively. 

Here is the rigorous definition of Pk for ^ £]3 • 2"^^', 4 • 2"'^[: 

f (1,0,0) if (yi,y2)G]0,2-2-'=[x]0,2-2-'=[ 



I3kir,yi,y2) ■■-- 



(-5, 0, 0) if (yi, y2) G]2 • 2-^ 4 • 2-'=[x]2 • 2-^ 4 • 2' 



(0,0,0) elsewhere on ]0, 4-2-'= [ X ]0, 4-2"^ [ 
Note that by construction 

div/3fc = on]0,4-2~''[x]0,4-2"^>]0,4-2-'=[. (4.13) 
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4.3.2. Construction of the vector field b. We now define the vector field b by 
using as a "building block" the vector field /3fc defined in § 4.3.1. We proceed 
in three steps. 

Step A: we extend /3fc to ]0, 2^~'^[xM^ by imposing that it is 2^~'^-periodic in 
both yi and y2, namely we set 

/3fc(r,yi+m22-^y2+n22-'=) := /3fc(r,yi,y2) (4.14) 

for every m,n £ Z and (2/1,2/2) G ]0,2^"^[ x ]0,2^"''[. We recall (4.13) and 
we observe that Pk is tangent at the surfaces yi = m2'^~^ and y2 = n2^~'^, 
m,n G Z. We therefore get 

div/3fc = on]0,22-^[xM2. (4.15) 

Step B: we define the vector field b{t,r,yi,y2) on the set A~ defined by (4.2). 
To this end, we introduce the decomposition 



]0,l[:=AAuU4, 



(4.16) 



fc=3 



where J\f is an ^^-negligible set and 

r ]V2,i 

4:= 



k fc-1 

i=3 i=3 



if A; 



if A; > 4. 



We then set 



bit,r,yi,y2) := /3fc(r + 1 - ^ 2^ ^yi,y2j 



)2-i 



in A , when r € I^. 



(4.17) 



i=3 



Some remarks are here in order. First, to illustrate the heuristic idea underlying 
definition (4.17) we focus on the behavior at time t = 1. The vector field 6(1, •) 
behaves like fS^ on the interval r G]1/2, 1[, like /S^ on the interval r G]1/4, 1/2[, 
like fS^ on the interval r G]1/8, 1/4[, and so on. In other words, as r — ^ 0^ the 
r-component of vector field b oscillates between the values 1, —5 and on a 
finer and finer "three-colors chessboard" . 

Second, we recall (4.15) and we observe that the vector field is continuous at 
the surfaces r = 1 — "^jZ^ 2^"-', k > 4. Hence, 



div 6 = on A 



(4.18) 



where the set A~ is the same as in (4.2). 

Step C: we define the vector field 6(t, r, yi, y2) on the set A+ defined by (4.1). 
To this end, we use the same decomposition of the unit interval as in (4.16) and 
we set 



Ki, r, yi, ys) := (1,0, 0)ln(t) + (-5, 0, 0)!^.^) in A+, when t £ 4. (4.19) 



(t) 
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In the previous expression, loit) denotes the characteristic function of the set 

k 
D(i):={(r,yi,y2):/3fc(t + l-5^22-i,yi,y2)-(l,0,0) = l} 

and l_B(j) is the characteristic function of the set 

k 
B{t) := {(r,yi,2/2) : /3fc(t + 1 - J^ 22-^,^1,2/2) • (1,0,0) = -s} 

i=3 

Hence, we have div 6 = on A+ since b is piecewise constant and tangent at the 
discontinuity surfaces. Next, we observe that by construction for any t G]0, 1[ 
the normal trace of b is continuous at the surface {(r, 2/1, 2/2) : r = t} C fi. By 
recalhng (4.18), we arrive at 

div6 = on]0, l[xJ]. (4.20) 

4.3.3. Proof of the regularity and trace properties. We first observe that, 
for every open and bounded set 0^= such that Q^, C il, we have b E 
L^ {{0,T[; BV {Q,^;W^)) since it is piecewise constant on il*. However, note that 
the BV regularity degenerates at the boundary r = 0. 

Next, we prove that Tr 6 = 1: we use the family of graphs {'^j-}re]o,i[ defined 
as in (4.7). By relying on Theorem 2.4, we infer that proving that Tr6 = 1 
amounts to show that there is a sequence r^ — t- 0"^ such that a^j. ^ 1 weakly* 
in L°°(]0, l[xD,.if^L]0, l[xL') as A; — )• +00, for every D open and bounded 
in R2. Hence, we can conclude by choosing r^ := 2^"'^. 

4.3.4. Construction of a nontrivial solution of the initial-boundary value prob- 
lem (1.4). To exhibit a nontrivial solution of (1.4) we proceed as follows: first, 
we give the rigorous definition, next we make some heuristic remark and finally 
we show that u is actually a distributional solution of (1.4). 

We recall the decomposition ]0, l[xQ = A"*" U A" U S, where A+, A~ and 
S are defined by (4.1), (4.2) and (4.3), respectively. We define the function 
u :]0, l[xil -^M by setting 

„(f ^ „ „ ^ ._ / 1 (t,r,yi,2/2) G A- and6(t,r,yi,j/2)-(l,0,0) = 1 . . 
^t,r,yi,y2).-| elsewhere in ]0,l[xO. ^^'^^^ 

The heuristic idea behind this definition is as follows. We have defined the vec- 
tor field b in such a way that, although b is overall outward pointing (namely, 
Tr6 > 0), there are actually countably many regions where b is inward point- 
ing (namely its r-component is strictly positive) which accumulate and mix at 
the domain boundary: these regions are represented by the dashed square in 
Figure 1. The function u is defined in such a way that u is transported along 
the characteristics (which are well-defined for a.e. (r, y) in the domain interior) 
and it is nonzero only on the regions where b is inward pointing. As a result, 
although b is overall outward pointing, it actually carries into the domain the 
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nontrivial function u. This behavior is made possible by the breakdown of the 
BV regularity of b at the domain boundary. 

We now show that n is a distributional solution of the initial-boundary value 
problem (1.4). First, we observe that u{t, •) ^ as t — )• 0"*" and hence the 
weak continuity of u with respect to the time implies that the initial datum is 
satisfied. 

We then set C := (n, bu) and we observe that Div C = on A^~ because C 
is identically there. Next, we observe that the vector field b is constant with 
respect to t in A~ and, by recalling that u is defined as in (4.21), we infer that 
u is also constant with respect to t in A~. Hence, showing that DivC = in 
A~ amounts to show that div (bu) = in A~. This can be done by relying on 
the same arguments we used to obtain (4.20). 

Finally, we observe that the normal vector to the surface S is (up to an 
arbitrary choice of the orientation) n := (l/v2, — l/v2, 0, 0). Hence, by con- 
struction the normal trace of C is zero on both sides of the surface S and hence 
DivCLcS = 0. This concludes the proof of Theorem 1.3. D 

4.4. Proof of Corollary 1.4. We first describe the heuristic idea underlying 
the construction of the vector field b. Loosely speaking, we proceed as in the 
proof of Theorem 1.3, but we modify the values of the "building block" /3fc 
on the subinterval r s]0, 2~'^[. Indeed, instead of defining /3fc as in Step 1 of 
§ 4.3.1, we introduce nontrivial components in the (yi,y2)-directions. These 
non-trivial components are reminiscent of the construction in Depauw [14] and 
the resulting vector field can be actually regarded as a localized version of 
Depauw's vector field. In particular, they enable us to construct a solution 
that oscillates between 1,-1 and and undergoes a finer and finer mixing as 

The technical argument is organized in two steps: in § 4.4.1 we introduce 
the "localized version" of Depauw vector field, while in § 4.4.2 we conclude the 
proof of Corollary 1.4. Before proceeding, we introduce the following notation: 

• Qk is the square (^1,2/2) G]0, 2-^[x]0, 2-*^[; 

• Sk is the square (2/1,2/2) e]0,2'^~^[x]0,2'^-''[. 

4.4.1. A localized version of Depauw [14] vector field. We construct the vector 
field afc, which is defined on the cell (r, 2/1,2/2) S]0, 2~'^[x(5fc. Also, for this 
construction we regard r as a time-like variable and we describe how a given 
initial datum evolves under the action of a^- The argument is divided into 
steps. 

Step 1: we construct the "building block" a^, which is defined on the square 
(2/1,2/2) G] - 2-2-'=,2-2-fc[x] _ 2-2-fc^2-2-'=[ by setting 

a.(yi,?/2)-| (22/2,0) I2/1KI2/2I. ^'-''^ 

Note that a^ takes values in M^, it is divergence free and it is tangent at the 
boundary of the square. 
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Figure 3. The action of the vector field a^ir, ■) on the solution 
Zk on the interval r s]0, 2~^~*^[. The solution Zj^ attains the 
values 1 and —1 on black and white regions, respectively. 




J/2 




Zk 



ak{r,-) 




Zk{2-'-'r) 



Step 2: we define the function Zj, : Q^ — )• M by considering the chessboard 
illustrated in Figure 3, left part. The function Zk attains the value —1 and 1 
on white and black squares, respectively. 

Step 3: we begin the construction of the vector field at :]0, 2~'^[x(5a; — > M^. 
If r g]0, 2~'^~^[, then ak{r, ■) is defined by setting 



a-kivi 



ak{r,yi,y2) 



2-^-^?/2 



-i-fc\ 



if (2/1,^2) G]2-2-^3 • 2-2-fc[x]2-2-^3 • 2 



-2-k( 



(4.23) 



(4.24) 



(0,0) elsewhere on Qfc- 

See Figure 3, central part, for a representation of the values attained by a^ on 
the interval r g]0,2"^~^'[. 

We term Zk the solution of the initial-boundary value problem 

drZk + divy{akZk) = in ]0, 2~'=[xQfc 
Zk = Zk at r = 0, 

where z^ is defined as in Step 2. Note that by construction diVyOfc = and 
therefore the first line of (4.24) is actually a transport equation. Hence, the 
value attained by the function z^ can be determined by the classical method 
of characteristics. In particular, the function Zfc(2~^~ , •) is represented in Fig- 
ure 3, right part, and it attains the values 1 and —1 on black and white squares, 
respectively. 

Step 4: if r g]2~^~''\ 3 • 2~^"'^[, then Qfc(r, •) is defined by setting 



afc(r,yi,2/2) = ak{yi 



^2^'"^y2 



J2- 



-2-/c^ 



if 



(2/1,2/2) G](i - 1)2-2-^ (i + 1) • 2-2-'^[x](i - 1)2-2-^ (j + 1) . i-^-f^i 



where i,j can be either 1 or 3. See Figure 4, central part, for a representation 
of the values attained by a^ on the interval r g]2~^~'^, 3 • 2^^~^[. Note that by 
construction divya^ = on ]0, 3-2~'^~^[xQk and hence the solution z^ of (4.24) 
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Figure 4. The action of the vector field afc(r, •) on the solution 
2fc on the interval r s]2~^~'^,3 • 2~'^~^[. The solution Zj^ attains 
the values 1 and —1 on black and white regions, respectively. 




Figure 5. The action of the vector field ak{r, •) on the solution 
Zk on the interval r €]3-2~^~'^, 4-2~^~'^[. The solution z^ attains 
the values 1 and —1 on black and white regions, respectively. 




evaluated at r = 3 • 2"^"*"' is as in Figure 4, right part: as usual, the black and 
white squares represent the regions where Zk{3 • 2~^~^,-) attain the values 1 
and —1, respectively. 

Step 5: if r g]3 • 2-2-'=, 4 • 2-2-'=[, then ak{r, ■) is again defined by (4.23). 
Hence, the values attained by Zk{2~^, •) are those represented in Figure 5, right 
part. 



4.4.2. Conclusion of the proof. Loosely speaking, the proof of Corollary 1.4 is 
concluded by combining the construction described in § 4.4.1 with the proof of 
Theorem 1.3. The argument is divided in four steps. 

Step A: we define the vector field /3fc and the solution Uk on {r,yi,y2) G 
]0,2-^[x5fc, where Sk is the square ]0, 22-^[x]0, 22-'=[. 



h{r,yi,y2) 
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We recall the definition of Bk provided by (4.9) and we set 

' (l,a,.(r,yi,y2)) (yi, ys) G]0, 2-^-[x]0, 2-'=[ 

(1, afc(r,yi - 2 • 2-^y2)) (yi, 2/2) G]2 • 2-^3 • 2-'=[x]0, 2" 

(1, ak{r, yi, 2/2 - 2 • 2-'=)) (yi, ys) G]0, 2-'=[x]2 • 2-^ 3 • 2" 
(l,afc(r,yi - 2 • 2-^y2 - 2 • 2~^)) {yi,y2) G]2-^3 • 2-^■[x]2-^3 • 

(-5,0,0) (2/1,2/2) GSfc 

(0, 0, 0) elsewhere in S^ 

Note that, basically, the definition of (3^ is obtained from (4.10) by changing 
the value of the vector field on D^ and inserting as a component in the (yi, y2)- 
directions the vector field a^ constructed in § 4.4.1. 
Also, we define the function u^ by setting 



2-^[ 



Ukir,yi,y2) ■-- 



Zk[r,yi,y2) 

Zfc(r,yl-2•2-^y2) 

^fc(r,yi,y2-2-2-'=) 



(yi,y2)G]o,2-'^[x]o,2-'=[ 

(yl,y2)G]2•2-^3•2-'=[x]0,2"^[ 
(yl,y2)G]0,2-'=[x]2•2-^3•2-'=[ 



Zk{r,yi - 2 • 2-^y2 - 2 • 2-^=) (yi,y2) G]2 • 2-^3 • 2-'=[x]2 • 2-^3 

elsewhere in S^, 

where z^ is the same function as in § 4.4.1. 

Step B: we define the vector field /J^ and the solution u^ for (r, yi,y2) € 
]2-^',22-^[x5fc. 

We set /3fc(r, yi, y2) := Pk{^, yi, y2), where /3fc denotes the same vector field as 
in § 4.3.1. The function Uk satisfies 

drUk + divyi^kUk) = 0. 

Since div^^^ = 0, the values attained by u^ for (r, yi,y2) g]2~^, 2^~'^[x5fc can 
be computed by the classical method of characteristics. To provide an heuristic 
intuition of the behavior of u^, we refer to Figure 1, center and right part, and 
we point out that u^ attains the value on white and black areas, while on 
dashed areas it attains the same values as in Figure 5, right part. 



•2-'=[ 



Step C: we extend /3k and Uk to ]0, 2' 



2-fcr 



by periodicity by proceeding 



as in (4.14). 

Step D: we finally define a vector field b and the function u. We recall the 
decomposition (4.16) and we define b as in (4.17) and (4.19), replacing /3fc with 
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/3fc. Also, we define u by setting 

, , _ J nfc(r yi,y2) in A", when r £ h 

By arguing as in the proof of Theorem 1.3, one can show that u and b satisfy 
requirements i), . . . , v) in the statement of Theorem 1.3 and that moreover 
Tr (bu) = 0. This concludes the proof of Corollary 1.4. D 
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